Time dilation of a bound half-fluxon pair 
in a long Josephson junction with a ferromagnetic insulator 
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The fluxon dynamics in a long Josephson junction with a ferromagnetic insulating layer is investi- 
gated. It is found that the Josephson phase obeys a double sine-Gordon equation involving a bound 
n fluxon solution, and the internal oscillations of the bound pair acting as a clock exhibit Lorentz 
reductions in their frequencies regarded as a relativistic effect in the time domain, i.e., time dilation. 
This is the complement to the Lorentz contraction of fluxons with no clock. A possible observation 
scheme is also discussed. 
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The Lorentz transformation in space-time plays a key 
role in the special theory of relativity, resulting in coun- 
terintuitive phenomena far from our everyday experience 
such as length contraction when an object is travelling 
close to the speed of light. The Lorentz contraction is 
difficult to confirm experimentally for energetic reasons 
related to macroscopic objects, i.e., macroscopic objects 
cannot be accelerated to a sufficiently high velocity for 
the Lorentz contraction to be observed. Only a few ex- 
amples have been demonstrated in specific systems using 
the Earth's rotation 0, 0, ■ elementary excitation 
with a quantum unit of magnetic fiux in long Josephson 
junctions, known as a fluxon, obeys sine-Gordon equa- 
tions involving the Lorentz invariance. The length con- 
traction of fiuxons in Josephson junctions was shown in 
experiments by recording the voltage pulse created by the 
fluxon motion 4] and by imaging the contraction of the 
fluxon- antifiuxon collision region in an annular Josephson 
junction with increasing fiuxon velocity by employing a 
low-temperature scanning electron microscope 5]. 

However, a single fiuxon cannot exhibit another rela- 
tivistic effect in the time domain, namely time dilation 
where moving clocks run slowly, because it does not have 
its own clock. As with any vibrating system, an internal 
degree of freedom in the system can serve as a moving 
clock counting time intervals. One candidate is a bound 
pair consisting of a fiuxon and an antifiuxon, and known 
as a breather. The breather moves along a Josephson 
transmission line with the internal oscillation in which 
the fluxon and the antifluxon are exchanging their posi- 
tions. This internal oscillation acts as a clock. However, 
the breather is inadequate for testing this effect because 
it is too fragile as regards external perturbations Q. In 
particular, the breather disappears even in a small dissi- 
pation because it is topologically equivalent to the vac- 
cum. This feature makes it extremely difficult to observe 
the effect of time dilation. Another possibility is the tt- 
TT kink pair solution in a double sine-Gordon equation 



(DSGE) in long isolated Josephson junctions [7]. This 
kink pair has a stable internal oscillation below the lower 
edge of the continuum phonon band 8]. Although this 
internal oscilation actually attenuates in the presence of 
dissipation, it is possible to excite the oscillation again 
on demand because the kink pair itself with a topological 
twist cannot decay into vacuum. However, it is difficult 
to observe tt-tt kink oscilations in isolated systems with 
a non-contact measurement. In this Letter, we propose 
a new type of Josephson transmission line that obeys a 
DSGE, namely a long Josephson junction with a ferro- 
magnetic insulator, in order to explore the time dilation 
of a bound pair of half fiuxons through their internal os- 
cillations. We also discuss an experimental scheme for 
observing time dilation in such a system. 

In the past decade, considerable attention has been di- 
rected towards Josephson tt junctions characterized by 
a minimum Josephson coupling energy at a phase dif- 
ference of TT, in relation to the physics of a Cooper pair 
with a finite momentum. The crossover between and 
TT junction states was demonstrated in superconductor- 
ferromagnet-superconductor (SFS) junctions as a func- 
tion of temperature [9i] and barrier thickness [lo| . In 
the vicinity of the crossover, it is conjectured that a 
second order component of Cooper-pair tunneling, i.e., 
sin(20) with (j> being the phase difference across the 
junction, becomes dominant in the current-phase rela- 
tion [H [H, 113, [13. This might significantly change 



the dynamics of the phase difference across the junction 
through the current-phase relation in a long Josephson 
junction. 

Consider two identical superconductors separated by a 
thin ferromagnetic insulator with thickness a lying in the 
xy plane as shown in Fig. 1. The dimension w in the y di- 
rection is assumed to be much smaller than the Josephson 
penetration depth, and thus the junction is considered to 
be one-dimensional [1]. Starting as usual from Maxwell 
equations that takes account of the spatiotemporal de- 
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FIG. 1: Schematic diagram of a long Josephson junction with 
a ferromagnetic insulator. 



pendence of the phase difference (j) across the junction 
due to electric field Ez and magnetic field By through 
the relations = {h/2ea)dt4' and By — {h/2ed)dx(f>, 
with dt = d/dt and dx = d/dx, we find 



{h/2e^ld)^l(i) = J + {he/2ea)dl(j) 



(1) 



where e is the dielectric constant of the ferromagnetic 
insulator. In contrast to a usual long Josephson junction, 
the current density in a long Josephson junction with a 
ferromagnetic insulator is now given as 



J — Jc{X sin (j) + sin 20) . 



(2) 



where Jc is the amplitude of the current density for the 
second harmonic and A is the ratio of the amplitude be- 
tween the first and second harmonics. When A ~ 0, Jc 
becomes the Josephson critical current density. Substi- 
tuting Eq. ([2|) into Eq. ([T]), we obtain 



di 



d^cj) = Asin0 



sin 2(f>, 



(3) 



where the coordinate x is normalized by v^Aj where 
Aj = y^h/4efidJc is the Josephson penetration depth for 
A = and the time t is normalized by ^/2ujJ^ where 
ujj = c/ Aj is the frequency of the Josephson plasma os- 
cillation with c = a/ a/defj,. The phase difference in an 
SFS hybrid junction thus obeys a DSGE, and constitutes 
one of the main results in this paper. 

A DSGE is not integrable and exact soliton solutions 
do not exist. However, it is possible to find a topologi- 
cally stable 27r-kink solution for a range of parameter A, 
i.e., A > 0. The potential energy density is expressed as 
Uj{(j)) = A(l — cos(f)) + i(I — cos 20). Here, the energy 
is normalized by Eq — 2^/2ejwXj where ej = ?iJc/4e is 
the specific Josephson coupling energy per unit area for 
A = 0. Since this potential has minima at = (mod 27r) 
for A > 0, there are 27r-kink solutions where the phase 
changes from to ±2it as x passes from — oo to oo. The 
shape of the 27r kink is determined by the competition 
between the potential energy density Uj{(j)) and the elas- 
tic energy density, C/e(0) — ^/2{dx(j))^, which measures 
the rigidity of the Josephson phase difference. 

For < A < 2, the potential energy density U,j also 
has an extra local minimum at = tt, in addition to the 
= minima. This minimum leads to a new situation 



where the 27r kink splits into two separate tt kinks ow- 
ing to the gain of the Josephson energy U,] against the 
elastic energy Ue- The equilibrium separation between 
two TT kinks is determined by the competition between 
the gain of C/j around = tt and the loss of Ue by the 
additional change in the phase difference. With a moving 
TT-TT kink pair, these compete in time, resulting in a new 
type of dynamics. The tt-tt kink pair exhibits an inter- 
nal oscillation related to the relative oscillations of two tt 
kinks around the equilibrium separation [isl . [l6l | . Figure 
[5] shows a typical profile of phase difference during the 
TT-TT kink oscillation. 




FIG. 2: (Color online). Phase profile during the tt-tt kink 
oscillation. Vertical (red) arrows denote the positions of tt 
kinks. Expansion and contraction of springs express the de- 
viation from equilibrium separation between tt kinks symbol- 
ically and holizontal (green) arrows denote the directions of 
restoring forces. 



The eigenfrequency of the tt-tt kink oscillation can be 
estimated by using the collective coordinate method [3] ■ 
We assume the solution to DSGE in the form 



2 tan 



.r(-+f) 



2 tan 



(4) 



where F = v2 + A is the slope of each tt kink and R 
denotes the separation between two tt kinks, which can 
be regarded as a collective coordinate that describes tt-tt 
kink oscillation. Substituting Eq. ((4]) into the Hamilto- 
nian of the DSGE, the potential energy and the inertial 
mass for the collective coordinate R are obtained as 



+ 2Ai?coth(ri?/2), 



;inh(rii) 



] coth2(ri?/2) 



(5) 
(6) 



The equilibrium separation Rq can be determined from 
the relation cosh^(Fi?o/2) = 1 + 2/A. Thus, the tt-tt 
kink oscillation can be regarded as the oscillation of a 
particle with a mass depending on its position around 
the equilibrium position Rq on the anharmonic potential 
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shown in Fig. [3l In a harmonic approximation we obtain 
the eigenfrequency wo(A) = \/Y" (Rq) /M(i?o) , where the 
prime denotes the differentiation with respect to R. This 
frequency of the tt-tt kink osciUation serves as a clock to 
measure the time. 
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FIG. 3: (Color online). Interaction potential between tt kinks. 

Since the DSGE is also Lorentz invariant, the period 
of the TT-TT kink oscillation should increase as the velocity 
of the pair increases as a consequence of the relativistic 
time delay 

T/To = {l-iv/cr}-'^\ (7) 

where T is the period of tt-tt kink oscillation at the finite 
velocity u, and Tg is the period at u = 0. However, in 
an actual junction, dissipation always exists and destroys 
Lorentz invariance. Scott studied the effects of dissipa- 
tion on the time dilation of a moving breather and showed 
that the frequency suffers Lorentz reduction, while the 
damping of oscillation is independent of velocity If 
this is also the case for tt-tt kink oscillation, a tt-tt kink 
pair can be used in testing the relativistic time dilation. 

To confirm the time dilation of a traveling kink pair, 
we performed a numerical simulation by using a standard 
difference approximation with the perturbed DSGE, 

dt(t) ~dl(f) + a dtcj) + \su\(j) + sin 2(j)^ h, (8) 

where a = I/lojRsC is a damping coefficient where Rg 
is the shunt resistance and C is the capacitance of the 
junction. is the direct bias current which keeps the 
velocity of the tt-tt kink pair. 

Figure |4] shows the numerical results for T/Tq as a 
function of the terminal velocity w of a pair, which is de- 
termined by the balance between damping (a) and driv- 
ing (/(,). Here we take A = 0.25 and assume that an 
initial bound pair is expressed by Eq. ^ with appropri- 
ate separations. The solid curve is simply the prediction 
obtained from the special theory of relativity, i.e., the rel- 
ativistic time delay, described by Eq. ([7]). From this, we 
can see that our numerical results for different strengths 
of dissipation are in good agreement with the theoreti- 
cal predictions. In addition, the values of Tq obtained 



from a numerical simulation are about 9.44 for a = 0.01 
and 9.46 for a = 0.05, which are close to the value of 
9.24 estimated from wo(0.25). For larger a, we have seen 
the agreement with the theory in the lower velocity re- 
gion where aT < 1, which means underdamping. These 
results indicate that the frequency of the tt-tt kink oscil- 
lation can be considered to be almost independent of a. 
Thus, the tt-tt kink oscillation is a promising candidate 
for use in testing the relativistic time dilation in super- 
conducting quantum circuits. 
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FIG. 4: Velocity dependence of the tt-tt kink oscillation period 
normalized by the period atv = for A = 0.25. The points are 
obtained by numerical simulations for different strengths of 
dissipation, while the solid curve is the theoretical expression 
of the time dilation, i.e., Eq. ([7]). 



Now let us discuss a possible scheme for observing the 
effect of the time dilatation of a tt-tt kink oscillation in 
a long SFS junction. In general, it is difficult to moni- 
tor the real-time dynamics of a tt-tt kink oscillation di- 
rectly, while it might be possible to determine the eigen- 
frequency of a propagating bound pair by using the reso- 
nance effect. The bound pair exhibits a forced oscillation 
under an external alternating current, leading to excita- 
tion of the TT-TT kink pair. In particular, resonance occurs 
when the frequency of the external current matches the 
frequency of the tt-tt kink oscillation. With resonance, 
the amplitude of the oscillation becomes large, and thus 
achieve a distinguishable value compared with the equi- 
librium width of the kink pair. It may be possible to 
observe the tt-tt kink pair by using an existing experimen- 
tal technique such as low-temperature scanning electron 
microscopy [5|, namely, by imaging the collision region 
between a TT-kink pair and an anti-TT-kink pair in an an- 
nular SFS Josephson junction. Thus, we can measure the 
frequency of the tt-tt kink oscillation by detecting the to- 
tal width of a tt-tt kink pair as a function of the frequency 
of the applied oscillating current. 

To study the resonance behavior, we performed a nu- 
merical simulation based on Eq. ([5]) including external al- 
ternating currents I^.^ sin cot together with the direct bias 
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current /b, where lex and lu are the normahzed ampH- 
tude and frequency, respectively. The parameters we use 
in the simulations are A = 0.25, /cx ~ 0.15 and a — 0.1. 




FIG. 5: (Color online). Amplitude of the tt-tv kink oscillation 
as a function of the period of external alternating current and 
the velocity of the bound fluxon pair. The (white) dotted line 
denotes the relation of the relativistic time dilation Eq. (O. 



Figure [5] shows the amplitude of the oscillation of the 
separation between tt kinks as a function of the normal- 
ized period of external alternating current 27r/wTo with 
To = 9.54 and the terminal velocity of a pair v/c. The 
amplitude is estimated from the half value of the dif- 
ference between the maximum and minimum separation 
of the kink pair, and is normalized by the width of a 
moving tt kink with velocity v, i.e., v^Aj^/l — {v/c)'^/r. 
We can easily find the resonance of the tt-tt kink oscil- 
lation as indicated by a (white) dotted line. This line 
shows the Lorentz dilation expressed by Eq. ([7|). In ad- 
dition, there is another resonance with 1/2 harmonic that 
has a period twice that of eigen oscillation as a result of 
a nonlinear effect, called nonlinear resonance [l3|. As 
seen in Fig. [3l a nonlinear interaction is induced be- 
tween vr kinks when the amplitude of the tt-tt kink os- 
cillation becomes large because of an increase in the ex- 
ternal field amplitude. In contrast to a linear interaction 
yielding a harmonic resonance, the nonlinear interaction 
gives rise to a subharmonic resonance at a double eigen- 
frequency. The 1/2-harmonic resonance corresponds to 
the cubic nonlinearity of the asymmetric potential shown 
in Fig. [3l There is also another resonance that is al- 
most independent of the pair velocity around a period 
of 0.5. This resonance corresponds to a homogeneous 
Josephson plasma oscillation whose period Tjp is deter- 
mined by the curvature of C/j around = minima, i.e., 
Tjp/To = 2TT/y/UjJfi)To = 27r/VXT2To ~ 0.44, where 
the prime denotes the differentiation with respect to cj). 

Finally, let us roughly estimate the parameters. In 
the vicinity of the 0-tt crossover, Josephson critical cur- 
rent would be quite small, say, two orders of magnitude 
smaller than that of usual junction [13, 14 1, and thus 



Aj '--'100 fim, wj ~10 GHz and a ~ 0.1. The amplitude 
of TT-vr kink oscillation at the resonance is of the order of 
10-100 iim, which is within the range of spatial resolution 
of the experiment Q . 

In summary, we have studied fluxon dynamics in a 
long Josephson junction with a ferromagnetic insulator. 
In such a system, the Josephson current-phase relation 
differs from that of a conventional junction, resulting 
in half-integer fluxons that obey a double sine-Gordon 
equation. The bound pair of the half-fluxons exhibits 
an internal oscillation that is unique to this system, and 
could be detected through this oscillation by using the 
resonance effect. This might provide evidence for the ex- 
istence of a second harmonic component in the current- 
phase relation of an SFS hybrid junction. In addition, we 
demonstrated numerically the time dilation of the bound 
half-fluxon pair by detecting the Lorentz reduction in its 
frequency as a function of pair velocity. This is the com- 
plement to the Lorentz contraction. Moreover, a quan- 
tized TT — TT kink pair might also provide an application 
for Josephson-based quantum computers as a new type 
of mobile qubit using Josephson tt states [l^ . 
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